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ABSTRACT 

We report on two programs to address differential evolution between the radio-loud 
and radio-quiet quasar populations at high (z > 4) redshift. Both programs entail 
studying the radio properties of optically-selected quasars. First, we have observed 
32 optically-selected, high-redshift (z > 4) quasars with the VLA at 6 cm (5 GHz). 
These sources comprise a statistically complete and well-understood sample. We detect 
t^- ' four quasars above our 3a limit of « 0.15 mJy, which is sufficiently sensitive to detect 

all radio-loud quasars at the probed redshift range. Second, we have correlated 134 
■ z > 4 quasars, comprising all such sources that we are aware of as of mid-1999, with 

FIRST and NVSS. These two recent 1.4 GHz VLA sky surveys reach 3a limits of 
approximately 0.6 mJy and 1.4 mJy respectively. We identify a total of 15 z > 4 quasars, 
of which six were not previously known to be radio-loud. The depth of these surveys 
Q-f does not reach the radio- loud/radio-quiet demarcation luminosity density (-L1.4 gh z = 

10 32 ' 5 h^Q ergss" 1 Hz -1 ) at the redshift range considered; this correlation therefore only 
provides a lower limit to the radio-loud fraction of quasars at high-redshift. The two 
programs together identify eight new radio-loud quasars at z > 4, a significant increase 
£> . over the seven currently in the published literature. We find no evidence for radio- 

•i-H . 

loud fraction depending on optical luminosity for —25 > Mb > —28 at z ~ 2, or for 
—26 > Mr > —28 at z > 4. Our results also show no evolution in the radio-loud 
fraction between z ~ 2 and z > 4 (-26 > M B > -28). 
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1. Introduction 

Already in the late 1960s it was becoming apparent that the rapid increase in quasar comov- 
ing space densities with redshift did not continue beyond z ~ 2. Several authors suggested the 
existence of a redshift cutoff beyond which a real decrease in quasar densities occurs, unrelated 
to observational selection effects (e.g.,Sandage 1972). Numerous studies have now shown that this 
is indeed the case. For example, Schmidt, Schneider, h, Gunn (1995a), applying the V/V max test 
to a well-defined sample of 90 quasars detected by their Lya emission in the Palomar Transit 
Grism Survey, conclude that the comoving quasar density decreases from z = 2.75 to z = 4.75. 
Warren, Hewett, & Osmer (1994) show that the quasar space density falls by a factor of a few 
between z = 3.3 and z = 4.0. Shaver et al. (1996b) show a similar decrease in space density exists 
for radio- loud quasars, implying that the decline is not simply due to obscuration by intervening 
galaxies. 

An unresolved question, however, is how the decline compares between the optically-selected 
and radio-selected populations. Two definitions are generally used to demarcate radio-quiet and 
radio-loud quasars. One criterion is the radio-optical ratio i? r _ of the specific fluxes at restframe 
6 cm (5 GHz) and 4400 A (Kellerman et al. 1989); radio-loud sources typically have R r ~ in the 
range 10 — 1000 while most radio-quiet sources have 0.1 < R r - < 1- However, Peacock, Miller, & 
Longair (1986) and Miller, Peacock, & Mead (1990) point out that R T - is physically meaningful 
as a discriminating parameter only if radio and optical luminosities are linearly correlated. This is 
apparently not the case; no linear correlation has been observed for radio-loud quasars (Stocke et al. 
1992) and the fraction of optically- faint, optically-selected quasars which are radio-loud would be 
higher if radio and optical properties were linearly correlated (Peacock et al. 1986; Miller et al. 
1990). 

The second definition, which can result in ambiguous taxonomy for some quasars, is to divide 
the populations at some restframe radio luminosity (Miller et al. 1990; Schneider et al. 1992). 
Different authors use slightly different luminosities to discriminate the populations; for what follows, 
we identify radio- loudness using the Gregg et al. (1996) cutoff value for the 1.4 GHz specific 
luminosity, L 1AGHz = 10 32 - 5 hV^ ergss" 1 Hz -1 ^ 10 24 WHz^sr 1 ). 

Differential evolution between the radio-quiet and radio-loud quasar populations would be a 
fascinating result, which could change (or challenge) our understanding of different types of active 
galactic nuclei (AGN). For example, Wilson & Colbert (1995) propose a model whereby radio-loud 
AGN are products of coalescing supermassive black holes in galaxy mergers; such a process would 
result in rapidly spinning black holes, capable of generating highly-collimated jets and powerful 
radio sources. Such a model could naturally explain a lag (if there is indeed one) between the 
appearance of powerful radio-loud quasars in the early Universe and the more common radio-quiet 
quasars. 

Several groups have reported on followup of optically-selected quasar surveys at radio frequen- 
cies (e.g., Kellerman et al. 1989; Miller et al. 1990; Schneider et al. 1992; Visnovsky et al. 1992; 
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Hooper et al. 1995; Schmidt et al. 1995b; Goldschmidt et al. 1999). Typically between 10% and 
40% of the quasars are detected in the radio. One notable exception is Kellerman et al. (1989) 
who, using the Very Large Array (VLA), detect more than 80% of 114 quasars observed from 
the Palomar-Green Bright Quasar Survey (BQS; Schmidt & Green 1983), a relatively low-redshift 
sample. Visnovsky et al. (1992) report on radio observations of 124 quasars in the redshift range 
1 < z < 3 selected from the Large Bright Quasar Survey (LBQS; Hewett, Foltz, & Chaffee 1995). 
The sample, chosen to complement the predominantly low-redshift BQS, leads them to conclude 
that the fraction of radio- loud quasars decreases with increasing redshift. A similar conclusion is 
reached by Schneider et al. (1992), who report on 5 GHz VLA observations of 22 optically-selected 
quasars at z > 3.1. La Franca et al. (1994), combining the results of several radio followup studies of 
optically-selected quasars, also concludes that the radio- loud fraction of optically-selected quasars 
decreases with increasing redshift; however, this conclusion is statistically robust only when the 
BQS sample is included. 

Hooper et al. (1995) report on 8.4 GHz VLA observations of an additional 132 LBQS quasars, 
supplementing the sample discussed by Visnovsky et al. (1992). Contrary to the results discussed 
above, Hooper et al. (1995) find that the radio- loud quasar fraction of the LBQS exhibits a peak 
around z ~ 1, is constant at ~ 10% until z ~ 2.5, and then increases substantially to z = 3.4, the 
highest redshift in the LBQS. They suggest the enhanced radio- loud fraction may be the result of 
increased radio-loud fraction at bright optical absolute magnitude (Mb Si —27.4). 

The most recent thorough analysis of the issue is presented by Goldschmidt et al. (1999), who 
report on 5 GHz VLA observations of 87 optically-selected quasars from the Edinburgh Quasar 
Survey. They combine these results with all published radio surveys of optically-selected quasar 
samples other than the BQS sample which has been shown to be incomplete by a factor of three 
(Goldschmidt et al. 1992). Miller, Rawlings, & Saunders (1993) suggest that this incompleteness is 
not random with respect to radio properties, implying that it is inappropriate to use the BQS for 
determinations of the radio- loud fraction. Goldschmidt et al. (1999) further fortify their sample 
by correlating several large quasar surveys with two recent, large-area, deep radio surveys: Faint 
Images of the Radio Sky at Twenty centimeters (FIRST; Becker, White, & Helfand 1995) and the 
NRAO VLA Sky Survey (NVSS; Condon et al. 1998). This work samples the optical luminosity- 
redshift plane more thoroughly than any previous analysis. Considering only quasars in the absolute 
magnitude interval —25.5 > Mb > —26.5, they find hints that the radio-loud fraction decreases 
with increasing redshift, though they note that the result is not statistically significant in the narrow 
Mb range considered. 

In short, the existence of differential evolution between the radio-quiet and radio-loud quasar 
populations remains uncertain. Some researchers find evidence that the radio-loud fraction of 
optically-selected quasars decreases with increasing redshift (e.g., Visnovsky et al. 1992), while 
others find no (or only marginal) evidence for evolution (e.g., La Franca et al. 1994; Goldschmidt 
et al. 1999), and yet other researches find that the radio- loud fraction increases with increasing 
redshift (e.g., Hooper et al. 1995). We attempt to re-address the issue by providing a more accurate 
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measure of the radio-loud fraction at z > 4. In §|2| we report on 5 GHz VLA observations of a 
sample of 32 z > 4 optically-selected quasars. The Schneider et al. (1992) 5 GHz VLA study 
includes 13 z > 4 quasars; we include these sources in our analysis. In §|3| we report on correlations 
of larger, more recent surveys of z > 4 quasars with deep, large-area radio surveys. §Q sets up the 
methodology for studying the radio-loud fraction, discussed in §|5|. 

Throughout we adopt the cosmology consistent with previous work in this field: an Einstein- 
de Sitter cosmology with Hq = 50 h^o kms" 1 Mpc -1 , Q = 1, and vanishing cosmological constant, 
A = 0. 

2. Observations 

2.1. VLA Sample Selection 

The 32 z > 4 quasars which comprise our VLA sample are listed in Table |]. The sample is 
from the Palomar multicolor survey (Kennefick, Djorgovski, & de Calvalho 1995; Djorgovski et al. 
1999) and represents all z > 4 quasars found from that survey at the time of the radio observations. 
Six of the quasars are equatorial objects found previously in a similar search by the Automated 
Plate Machine (APM) group (Irwin, McMahon, & Hazard 1991) which fall within the Palomar 
multicolor survey selection criteria. These 32 sources comprise a statistically complete and well- 
understood sample of optically-selected z > 4 quasars; they should be an unbiased set as far as the 
radio properties are concerned. 

Calculation of the absolute B magnitudes listed in Table |] is discussed in §§|. 

2.2. VLA Observations 

The VLA observations were made on UT 1997 March 14, from 01:30 to 17:30 LST. The targets 
were observed at 4835 and 4885 MHz. Each frequency was observed with 50 MHz bandwidth and 
both polarizations. The configuration was 'B', with a maximum baseline of approximately 10 km 
and a resolution of ~ l'.'3. Observations were made in nodding mode, in which we oscillated between 
40 s observations of a calibrator and 240 s observations on a target. This procedure allowed us 
to better track and remove atmospheric fluctuations. Calibration and imaging followed standard 
procedures. Each target was observed for 30 — 45 min, including calibration. The rms noise was 
typically between 30 and 60 /iJy 

We detect four quasars from this sample in the VLA imaging, where a match is liberally defined 
to correspond to a radio source lying within 10" of the optical source. In fact, the four matches 
all have optical-to-radio positional differences less than 2", implying that the identifications are 
unlikely to be spurious. Three sources represent newly identified high-redshift, radio- loud quasars. 



- 5 - 



Table 1: Optical and Radio Properties of the High-Redshift VLA Sample 







r 


R.A. 


Dec. 


Mb 


SsGHz 


Quasar 


z 


(mag) 


(J2000) 


(J2000) 


(mag) 


(mJy) 



PSS0003+2730 


4 


.26 


19. 





00 


03 


23 


.10 


+27 


30 


23. 


.0 


-27, 


54 


< 


0.171 


PSS0030+1702 


4 


.28 


19. 


3 


00 


30 


16 


.39 


+17 


02 


40. 


.2 


-27, 


27 


< 


0.141 


PSS0034+1639 


4 


.38 


19. 


.5 


00 


34 


54. 


.82 


+16 


39 


19. 


.5 


-27, 


23 


< 


0.171 


PSS0059+0003 


4 


.16 


19. 


.5 


00 


59 


22 


.80 


+00 


03 


01. 


.2 


-26, 


90 


< 


0.138 


PSS0106+2601 


4 


.32 


19. 


.4 


01 


06 


00. 


.80 


+26 


01 


02. 


.5 


-27, 


23 


< 


0.144 


BRI0103+0032 


4 


.43 


18. 


6 


01 


06 


19 


.14 


+00 


48 


23. 


.1 


-28, 


22 


< 


0.216 


PSS0117+1552 


4 


,24 


18. 


6 


01 


17 


31. 


.17 


+15 


52 


16. 


3 


-27, 


91 


< 


0.183 


PSS0132+1341 


4 


.16 


19. 


.3 


01 


32 


09 


.85 


+13 


41 


39. 


.5 


-27, 


,10 


< 


0.165 


PSS0134+3307 


4 


.52 


18. 


.8 


01 


34 


21. 


.60 


+33 


07 


56. 


.5 


-28, 


18 


< 


0.165 


PSS0137+2837 


4 


.30 


19. 





01 


37 


12. 


.30 


+28 


37 


34. 


.5 


-27, 


60 


< 


0.171 


PSS0152+0735 


4 


.07 


19. 


6 


01 


52 


11. 


.10 


+07 


35 


50. 


.2 


-26, 


70 


< 


0.165 


BRI0151-0025 


4 


.20 


18. 


.9 


01 


53 


39 


.61 


-00 


11 


04. 


9 


-27, 


55 


7.96 


±0.110 


BRI0241-0146 


4 


.04 


18. 


.2 


02 


44 


01 


,85 


-01 


34 


03. 


.0 


-28, 


07 


< 


0.162 


PSS0244-0108 


4 


.00 


19. 


.0 


02 


44 


57 


.18 


-01 


08 


08. 


.6 


-27, 


25 


< 


0.180 


PSS0248+1802 


4 


.43 


18. 


4 


02 


48 


54. 


.28 


+18 


02 


49. 


3 


-28, 


,42 


0.39 


± 0.052 


PSS0747+4434 


4 


.42 


18. 


1 


07 


47 


50. 


,00 


+44 


34 


16. 


.0 


-28, 


71 


< 


0.126 


BRI0952-0115 


4 


.43 


18. 


7 


09 


55 


00 


,08 


-01 


30 


06. 


.9 


-28, 


,12 


< 


0.120 


BRI1013+0035 


4 


.38 


18. 


8 


10 


15 


48. 


.96 


+00 


20 


19. 





-27, 


93 


< 


0.114 


PSS1048+4407 


4 


.45 


19. 


2 


10 


48 


46 


62 


+44 


07 


12. 


7 


-27, 


.66 


< 


0.159 


BRI1050-0000 


4 


.29 


18. 


8 


10 


53 


20. 


.43 


-00 


16 


49. 


6 


-27, 


78 


10.2 


± 0.044 


PSS1057+4555 


4 


.12 


17. 


7 


10 


57 


56. 


39 


+45 


55 


51. 


9 


-28, 


65 


< 


0.102 


PSS1159+1337 


4 


.08 


18. 


5 


11 


59 


06 


.48 


+13 


37 


37. 


.7 


-27, 


.81 


< 


0.123 


PSS1317+3531 


4 


.36 


19. 


1 


13 


17 


43. 


.21 


+35 


31 


31. 


1 


-27, 


60 


< 


0.132 


PSS1430+2828 


4 


.30 


19. 


3 


14 


30 


31 


.63 


+28 


28 


34. 


3 


-27, 


30 


< 


0.108 


PSS1435+3057 


4 


.35 


19. 


3 


14 


35 


23. 


52 


+30 


57 


22. 


5 


-27, 


.38 


< 


0.108 


PSS1438+2538 


4 


.25 


19. 


5 


14 


38 


35. 


.58 


+25 


38 


32. 


1 


-27, 


.02 


< 


0.105 


PSS1443+2724 


4 


.42 


19. 


3 


14 


43 


31. 


.17 


+27 


24 


37. 





-27, 


.51 


< 


0.114 


PSS 1456+2007 


4 


.25 


19. 


5 


14 


56 


28. 


.97 


+20 


07 


27. 


.1 


-27, 


.02 


1.30 


± 0.049 


PSS1618+4125 


4 


.21 


19. 


6 


16 


18 


22. 


73 


+41 


25 


59. 


7 


-26, 


86 


< 


0.129 


PSS1646+5514 


4 


.05 


18. 


1 


16 


46 


56 


.29 


+55 


14 


46. 


7 


-28, 


.18 


< 


0.105 


PSS1721+3256 


4 


.03 


19. 


2 


17 


21 


06. 


.74 


+32 


56 


35. 


7 


-27. 


06 


< 


0.138 


PSS2122-0014 


4. 


.18 


19. 


1 


21 


22 


07. 


.50 


-00 


14 


45. 





-27, 


32 


< 


0.243 



Notes. — Radio positions are listed for detected sources; all other positions are from the optical. Where no radio 
source was detected within 10" of the quasar optical position, we quote the 3cr noise measured at the center of the 
VLA map as an upper limit to the 5 GHz flux density. See text for description of Mb determination. PSS quasars are 
from the Palomar multicolor survey. BRI quasars are from the APM survey. Optical coordinates for BRI 1050— 0000 
are from Storrie-Lombardi et al. (1996), which are « 50" offset from the coordinates reported by Smith et al. (1994). 

The quasar BRI1050— 0000 was previously identified as a 10.6 ± 0.2 mJy source at 4.8 GHz by 
McMahon et al. (1994). The source PSS1618+4125 has a 4a radio detection 16'.'0 away from the 
optical identification; we do not consider this a positive match. 
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Table 2: Optically-Selected, High-Redshift Quasars Detected by FIRST 







r 




R.A. 


Dec. 


M B 


Sl.4GHz 


Offset 




Quasar 


z 


(mag) 




(J2000) 


(J2000) 


(mag) 


(mJy) 


(arcsec) 


Ref 


BRI0151-0025 


4.20 


18.9 


01 


53 39.616 


-00 11 04.63 


-27.55 


4.75 + 0.15 


0.90 


1-* 


PSS1057+4555 


4.12 


17.7 


10 


57 56.324 


+45 55 53.04 


-28.65 


1.38 + 0.13 


2.13 


2-* 


FIRST J141045+34009 


4.35 


19.6 


14 


10 45.758 


+34 09 09.68 


-27.10 


2.07 + 0.15 


0.84 


3 


GB1428+4217 


4.72 


20.9 


14 


30 23.740 


+42 04 36.52 


-26.44 


215.62 + 0.15 


0.41 


4 


FIRST J 145224+335424 


4.12 


20.4 


14 


52 24.254 


+33 54 24.64 


-25.95 


6.95 + 0.13 


3.02 


3 


FIRST J171356+421808 


4.23 


19.0 


17 


13 56.150 


+42 18 08.64 


-27.49 


2.80 + 0.15 


0.00 


3 



Notes. — Newly identified high-redshift, radio-loud quasars are indicated with a *. Positions are from the FIRST 
catalog. References: (1) Smith et al. (1994); (2) Kennefick et al. (1995); (3) Stern et al. (1999); (4) Hook & McMahon 
(1998). 

3. Optically-Selected z > 4 Quasars in FIRST and NVSS 

To augment the targeted VLA sample discussed above, we have correlated an updated list 
of all 134 z > 4 quasars known to us in mid-1999 with two recent, large-area radio surveys. The 
FIRST survey (Becker et al. 1995), which is still in production mode, is a radio map of the northern 
celestial sphere at 21 cm (1.4 GHz), reaching a typical limiting flux density of 1.0 mJy (5a). The 
1.4 GHz NVSS survey (Condon et al. 1998) covers a larger area to a shallower flux density limit; 
the 5<T limit of NVSS is 2.25 mJy. For luminosity distance and radio spectral index a (S u oc v a ), 
the restframe specific luminosity L v = Aird 2 L S v / (\ + z) 1+a , where both L v and S v are measured 
at the same frequency. At z = 4 and for a typical quasar spectral index of a = —0.5, the FIRST 
survey reaches a 3d limiting 1.4 GHz specific luminosity of log Li. 4 (h^ ergss -1 Hz" 1 ) = 32.6. The 
comparable limit for the NVSS survey is logLi.4 (/i^q 2 ergss -1 Hz" 1 ) = 32.9. Therefore, using the 
radio luminosity definition of radio loudness, these surveys incompletely sample radio- loud quasars 
at high redshift. 

Of the 134 z > 4 quasars, 51 reside in portions of the celestial sphere observed thus far by 
FIRST. We deemed a radio detection in FIRST to be associated with a high-redshift quasar if it 
lay within 10" of the quasar optical position. Six high-redshift, radio-loud quasars were identified, 
of which four were previously known (see Table |). One, GB1428+4217 (Hook & McMahon 1998), 
is the highest-redshift radio-loud quasar known currently (z = 4.72); it was initially identified on 
the basis of being a strong, flat-spectrum radio source coincident with an unresolved, red, stellar, 
optical source. The three sources with VLA nomenclature were initially identified by correlating 
the FIRST survey with the digitized 2nd generation Palomar sky survey (Stern et al. 1999). 
Two new, high-redshift, radio-loud quasars have been identified: BRI0151— 0025 at z = 4.20 and 
PSS1057+3409 at z = 4.12. The former was also identified in the targeted 5 GHz survey discussed 
above (§|j), implying a radio spectral index a = 0.41. The latter was undetected in the targeted 
5 GHz survey, implying an unusually steep spectral index, a < —2.05 (3 a) or significant variability 
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Table 3: Optically-Selected, High-Redshift Quasars Detected by NVSS 







y 


R.A. 


Dec. 


Ma 


5 

.4GHz 


Offset 




Oil fltiflT" 






( T2000") 


( 12000") 








Ref 


pnoo974-0595 


4.10 


21.5 


00 29 50 18 

UU 0\J . io 


+05 42 1 2 9 


—24.82 


4.8 ± 0.5 


8.47 


1— * 


PSS0121+0347 


4.13 


17.9 


01 21 26.22 


+03 47 06.3 


-28.46 


78.6 ± 2.4 


0.76 


2-* 


BRI0151-0025 


4.20 


18.9 


01 53 39.71 


-00 11 12.1 


-27.55 


3.0 + 0.5 


7.57 


3-* 


RXJ1028.6-0844 


4.28 


18.9 


10 28 38.82 


-08 44 38.6 


-27.67 


269.8 + 8.1 


22.4 


4 


BRI1050-0000 


4.29 


18.8 


10 53 20.28 


-00 16 51.0 


-27.78 


9.7 + 0.6 


3.04 


5.6 


BRI1302-1404 


4.04 


18.6 


13 05 25.16 


-14 20 41.0 


-27.67 


20.6 + 0.8 


0.61 


7-* 


FIRST J141045+340909 


4.36 


19.6 


14 10 45.87 


+34 09 22.4 


-27.10 


3.5 + 0.15 


13.5 


8 


GB1428+4217 


4.72 


20.9 


14 30 23.69 


+42 04 36.8 


-26.44 


211.1 + 6.3 


1.10 


9 


FIRST J 145224+335424 


4.12 


20.4 


14 52 24.11 


+33 54 23.0 


-25.95 


8.6 + 0.5 


2.42 


8 


GB1508+5714 


4.30 


18.9 


15 10 02.90 


+57 02 43.3 


-27.70 


202.4 + 6.1 


22.3 


10 


FIRST J171356+421808 


4.23 


19.0 


17 13 56.37 


+42 18 06.3 


-27.49 


2.9 + 0.5 


3.40 


8 


GB1713+2148 


4.01 


21.0 


17 15 21.16 


+21 45 34.2 


-25.25 


447.2 ± 15.0 


2.02 


9 


RXJ1759.4+6638 


4.32 


19.0 


17 59 26.82 


+66 38 43.5 


-27.63 


3.6 + 0.7 


11.6 


11-* 


PC2331+0216 


4.10 


20.0 


23 34 32.1 


+02 33 16.7 


-26.32 


2.7 + 0.5 


6.51 


12,13 



Notes. — Newly identified high-redshift, radio-loud quasars are indicated with a *. Positions are from the NVSS 
catalog. References: (1) Schneider et al. (1997); (2) Djorgovski et al. (1999); (3) Smith et al. (1994); (4) Zickgraf 
et al. (1997); (5) Smith et al. (1994); (6) McMahon et al. (1994); (7) Storrie-Lombardi et al. (1996); (8) Stern et 
al. (1999); (9) Hook & McMahon (1998); (10) Hook et al. (1995); (11) Henry et al. (1994); (12) Schneider, Schmidt, 
& Gunn (1989); (13) Schneider et al. (1992). Optical coordinates for BRI1050-0000 are from Storrie-Lombardi et 
al. (1996), which are w 50" offset from the coordinates reported by Smith et al. (1994). 

on the time scale corresponding to the epochs of the two radio surveys. 

We have also correlated the list of z > 4 quasars with the NVSS catalog. 129 of the 134 quasars 
are in portions of the sky surveyed by NVSS. A radio detection in NVSS was deemed associated 
with a high-redshift quasar if it lay within 30" of the quasar optical position. We found a total of 
14 radio identifications (see Table 3); 6 of these quasars were initially discovered on the basis of 
their radio emission (Hook et al. 1995; Hook & McMahon 1998; Stern et al. 1999) and three were 
previously known to be radio-loud (Schneider et al. 1992; McMahon et al. 1994; Zickgraf et al. 
1997). Five new high-redshift, radio- loud quasars have been identified, of which two were identified 
in the targeted 5 GHz survey (§g): BRI0150-0025 (a = 0.77) and BRI1050-0000 (a = 0.04). 

Although 30" is a rather large matching radius, the surface density of sources in NVSS is 
~ 60 deg -2 , so that the probability of a chance coincidence of an NVSS source with an optical 
source is only 0.013. We therefore expect that 1.7 of our radio identifications are spurious. We 
note that of the four radio-optical offsets larger than 10", the three largest correspond to known 
radio-loud quasars, illustrating that the 30" search radius was appropriate. The fourth largest 
radio-optical offset belongs to RXJ1759. 4+6638, an X-ray selected quasar identified from deep 
ROSAT observations near the North Ecliptic Pole (Henry et al. 1994). NVSS reports a 3.6 mJy 
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source 11'' 6 away. Sensitive VLA 1.5 GHz C-array mapping of this region by Kollgaard et al. 
(1994) finds no object brighter than 0.5 mJy associated with the source, apparently at odds with 
the NVSS results. 

We are aware of only one z > 4 radio-loud quasar not included by these exercises: PKS1251— 407 
at z = 4.46 (Shaver, Wall, & Kellerman 1996a) which resides too far south to be within the VLA 
surveys. 

The list of z > 4 quasars comes from diverse sources, limiting the significance of these detections 
to studies of the quasar luminosity function at high redshift. However, 107 of the 129 z > 4 quasars 
were initially identified from optical surveys, either the Palomar multicolor survey (Djorgovski et al. 
1999), the APM survey (Irwin et al. 1991), or the long-term program of Schneider et al. (1989, 
1997). These surveys should be largely unbiased with respect to radio properties. In the following 
analysis, FIRST/NVSS detections of these sources will provide lower limits to the radio- loud quasar 
fraction at high-redshift. 

4. Results 

4.1. 1.4 GHz Specific Luminosity 

The specific luminosity at 1.4 GHz is straight-forward to calculate. For S u oc v a at radio 
frequencies, 

L u (, 1 )=4ndl {i + 1 z)1+a {^JSM) (1) 

where L v {y\) is the specific luminosity at rest-frame frequency v\ and S u {v2) is the flux density at 
observed frequency 1/2 • For the adopted cosmology, chosen to be consistent with previous work in 
this field, the luminosity distance 

d L = (2c/H Q )(l + z-Vl + ~z). (2) 

In Fig. H] we present the results of the radio surveys, plotting L1.4 GHz against redshift, with the 
radio luminosity cutoff between radio-loud and radio-quiet quasars indicated by a horizontal dashed 
line. Again, consistent with previous work in this field, we take a = —0.5; this is a typical value 
for radio- loud quasars. 

4.2. Absolute B Magnitude 

In order to compare our results at z > 4 with previous results, we require transforming the 
apparent r magnitude, m r (A Q b s ~ 6400 A), to a restframe absolute B magnitude, Mb (A rest ~ 4400 
A). We follow the methodology of Kennefick et al. (1995), in which a model quasar spectrum is 
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Fig. 1. — Li.4 GHz plotted against z for z > 4 quasars. Legend relates symbol to radio source 
catalog. 3o" limits of NVSS and FIRST surveys are illustrated with dotted lines. PSS/BRI refers 
to the 5 GHz deep imaging discussed in §|2|. SSG92 refers to the Schneider et al. (1992) sample. 
Non-detections are indicated with open symbols. Horizontal dashed line indicates adopted dividing 
line between radio-loud quasars and radio-quiet quasars. The k-corrections have been calculated 
assuming a = —0.5. 



used to relate m r to the apparent magnitude at wavelength 1450(1 + z) A. A power-law optical 
spectrum is then assumed longward of Lya to relate this to Mb- 



The absolute magnitude M is given by (e.g., Hogg 1999), 

Mab(A) = m AB (A) - 51og(d L /10 pc) + 2.5 log [(1 + z) 



L v [X/{l + z) 



L v {\) 



(3) 



where both the apparent magnitude m and the absolute magnitude M are measured at the same 
wavelength. The middle term on the right is called the distance modulus and the last term is 
the k-correction. Magnitudes are referred to the AB system (Oke 1974), defined as mAB(A) = 
—2.5 log S u (\) — 48.60, where apparent flux density S u is measured in units of ergscm^ 2 Hz -1 . 
Absolute magnitudes at differing wavelengths are related by 



Mab(Ai) - M AB (A 2 ) = -2.5 log 



MAi) 



(4) 
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Setting Ai = 4400 A and A2 = 1450(1 + z) A, the above equations are combined to yield 

L„(1450) 



M AB (4400) = ?7tAB [1450(1 + z)} -5log(d L /W pc) + 2.5 log 



1 + * 



^,(4400) 



(5) 



We assume a power-law optical spectrum longward of Lya, L u oc i/ a °v*. For a standard quasar op- 
tical spectral index a op t = —0.5 (e.g.,Richstone & Schmidt 1980; Schneider et al. 1992), consistent 
with previous work in this field, the offset between Vega-based Mb and ^4-B-system Mab(4400) is 
given by Kennefick et al. (1995), 

M B = M AB (4400) + 0.12. (6) 

We note that the adopted optical spectral index is likely too flat. Preliminary results from near- 
infrared imaging of z > 4 quasars to target the rest-frame B show substantial scatter in spectral 
slopes, but with an average slope slightly steeper than the standard (lower-redshift) value (Djor- 
govski et al. , in preparation). The final expression relating Mb to 77j A b[1450(1 + z)] is then 

M B = m AB [1450(l + z)] - 51og(d L /10 pc) + 2.51og(l + z) + 1.21a opt + 0.12. (7) 



Our observable is the (Vega-based) r magnitude, m r . An analytic expression relating the 
apparent magnitude at 1450(1 + z) A to m r is not possible for the redshift range we are considering 
as Lya emission and the (redshift-dependent) Lya forest traverse the r-band for 4 < z < 5. Instead 
we compute this relation, or k-correction, 

A; = m AB [1450(l+z)]-m AB (r) (8) 

using the Francis et al. (1991) composite (LBQS) quasar spectrum modified by the Madau (1995) 
model of the hydrogen opacity of the Universe. We restrict this calculation to the Lya forest 
(ignoring the other Lyman series absorptions). For A Q b s < 1216(1 + z) A, Madau (1995) finds the 
optical depth of the intergalactic medium is well-represented as the sum of Lya-forest absorptions 
and hydrogen absorption associated with metal line systems: respectively, 

/ X \ 3 - 46 / X \ 168 

r cfr = 0.0036 ( ^L) + 0.0017(^1 . (9) 
V 1216/ V 1216/ 

This model well describes the observed flux decrements at the redshifts considered. 

We assume that the Lya forest is sparse in the LBQS composite, relative to the thick jungle we 
are calculating at z > 4 and scale the LBQS composite by e~ TcS shortward of Lya. The k-correction 
is calculated by convolving the resultant, redshift-dependent model spectra with the Gunn r filter 
(see Fig. ||). We show the same curve for the Kron R filter as well; these k-corrections match 
reasonably well with the calculations of Kennefick, Djorgovski, & Meylan (1996) for a different 
model quasar spectrum and different parameterization of the Lya forest. At the desired accuracy 
of these calculations, it is appropriate to adopt the flat-spectrum (in S u ) relation between AB and 
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Fig. 2. — The k-correction between apparent AB magnitude at 1450 (1 + z) A and apparent 
AB magnitude through Gunn r filter, mp&(r), plotted as a function of redshift (solid line). The 
similar relation is also plotted for the Kron R filter (dashed line). As described in the text, the k- 
correction was calculated for the Francis et al. (1991) composite LBQS quasar spectrum modified 
by the redshift-dependent strength of the Lya forest, as determined from the model of Madau 
(1995). 



Vega-based r magnitude, r = m r « jjiab (?") — 0.21, so that our final expression relating Mb to the 
observed m r is 

M B = m r - 5 log(d L /10 pc) + 2.5 log(l + z) + 1.21a opt + k cS , (10) 
where k c g = k + 0.33. 

This method of determining Mb is slightly different from the approach used by Schneider et al. 
(1992) in their VLA followup of 22 optically-selected quasars at z > 3.1. Schneider et al. (1992) use 
the extinction-corrected spectrophotometric flux density at 1450 x (1 + z) A and assume an optical 
spectral energy index of —0.5 to calculate the absolute B magnitude. We have recalculated Mb 
for their sample using our method. The root-mean-square (rms) difference between the methods is 
0.30 mag for the 13 quasars above z = 4 in their sample. If we only consider the nine quasars at 
4.0 < z < 4.4, the rms difference between the calculated absolute B-bwd magnitudes is 0.23 mag. 
This is comparable to the photometric uncertainties. In what follows, we use the revised Mb for 
the Schneider et al. (1992) sample. Absolute B magnitudes should be correct to ~ 0.3 magnitudes. 
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Fig. 3. — Location of the z > 4 quasars in the optical luminosity-redshift plane. PSS/BRI refers 
to the 5 GHz deep imaging discussed in §g. SSG92 refers to the Schneider et al. (1992) sample. 
Non-detections are indicated with open symbols. The z > 4 FIRST and NVSS non-detections are 
not presented. 



4.3. Coverage of the Optical Luminosity Redshift Plane 

Fig. [H summarizes the location of the z > 4 quasars in the optical luminosity-redshift plane. 
Consonant with their extreme distances, these quasars have extremely bright absolute magnitudes, 
Mb ^ —26. To probe evolution in the radio-loud quasar fraction, we require comparison with a 
sample of comparably luminous quasars at lower redshift. Fig. || illustrates such a sample: we have 
augmented the z > 4 quasars with several lower-redshift surveys from the literature, summarized 
below. As before, k-corrections have been calculated assuming that both the radio spectral index a 
and the optical spectral index a opt are equal to —0.5. We do not include radio surveys of optically- 
selected quasars which have few objects at our target optical luminosity (e.g., the Edinburgh quasar 
sample considered in Goldschmidt et al. 1999). 

Miller sample: Miller et al. (1990) report on sensitive 5 GHz VLA observations of 105 quasars 
at 1.8 < z < 2.5. They provide apparent magnitudes at the restframe wavelength of 1475 A. 
Following La Franca et al. (1994) and Goldschmidt et al. (1999), we transform to B using uib = 
mi475 + 0.23; this is correct at the average redshift of the sample for the assumed optical spectral 
index. The typical limiting flux density of this survey is Si.4GHz = 0.5 mJy (3a). For z ~ 2, 
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M B (mag) 

Fig. 4. — Location of all quasars considered here in the optical luminosity-redshift plane. PSS/BRI 
refers to the 5 GHz deep imaging discussed in §§. SSG92 refers to the Schneider et al. (1992) 
sample. MPM90 refers to Miller et al. (1990). LaF94 refers to La Franca et al. (1994). LBQS 
refers to the 8.4 GHz VLA imaging of LBQS quasars reported in Visnovsky et al. (1992) and 
Hooper et al. (1995). Non-detections are indicated with open symbols. The z > 4 FIRST and 
NVSS non-detections are not presented. Vertical dotted lines indicate the optical luminosity range 
considered in our analysis. 

this translates to a radio luminosity of Li.4GHz — 1-5 x 10 32 H§q ergss _1 Hz _1 , well below the 
radio-loud cutoff. 

La Franca sample: La Franca et al. (1994) report on 5 GHz VLA observations of 23 optically- 
selected quasars with B < 19.4 at intermediate redshift (0.8 < z < 3.4). The average 3a limits to 
detections are ~ 0.12 mJy, corresponding to Li.4gh z — 8 x 10 31 /i^~ 2 ergss -1 Hz _1 at z ~ 3. This 
is well below the radio-loud cutoff. We use the absolute luminosities provided in the paper. 

LBQS sample: 256 quasars from the LBQS have been observed at 8.4 GHz by the VLA by 
Visnovsky et al. (1992) and Hooper et al. (1995). This large survey, sampling 0.2 < z < 3.5 and 
—22.5 > Mb ^ —29, suffers from strongly correlated luminosity and redshift, as expected for a 
flux-limited sample: at a given redshift, Mb only spans « 1.5 magnitudes. We use the absolute 
magnitudes provided in the papers. The median 3a noise limit of 58.4GHz = 0.29 mJy corresponds 
to Li.4GHz — 2.5 x 10 32 h$Q ergss _1 Hz _1 , again below the radio-loud cutoff. We do not include 
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Fig. 5. — Location of all quasars considered here in the optical luminosity-radio luminosity plane. 
Symbols are as in Fig. ||; again, non-detections are indicated with open symbols. Horizontal dashed 
line refers to the cutoff between radio- loud and radio-quiet quasars. A few sources with elevated 
noise levels are classified as radio-loud. Vertical dotted lines indicate the optical luminosity range 
considered in our analysis. 



the additional 103 LBQS quasars with 8.4 GHz VLA observations discussed in Hooper et al. (1996) 
as that sample is largely optically less-luminous, lower-redshift quasars. 

Fig. [B| illustrates the location of the quasars in Mb — £i.4GHz space. For the purposes of 
the following analysis, we assume non-detections have radio fluxes equal to their 3a noise value. 
We omit the new z > 4 radio-loud quasars detected by FIRST and NVSS in the following anal- 
ysis, as those surveys are insufficiently sensitive to reach the radio-loud cutoff at z > 4. The 
detection limits of the other surveys are sufficiently deep that few non-detections are above the 
radio-loud/radio-quiet boundary; we conservatively classify those few sources as radio- loud be- 
low. Our total sample is 428 quasars spanning 0.2 < z < 4.7, —22.7 > Mb > —28.7, and 
30.08 < logL^GHz^ergss^Hz- 1 ) < 35.7. 



5. Radio-Loud Fraction 
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Fig. 6. — Fraction of radio-loud quasars as function of absolute magnitude. Two redshift ranges 
are plotted with 1 mag wide bins. No evidence for luminosity dependence in the radio-loud fraction 
is evident from this limited data set. 



5.1. Radio-Loud Fraction as a Function of Optical Luminosity 

We first consider if we detect optical luminosity dependence in the radio-loud fraction. We con- 
sider two redshift ranges where we have large samples of quasars, and divide the sample into absolute 
magnitude bins where the quasars are approximately smoothly distributed. Fig. || shows this frac- 
tion for each redshift bin. Error bars shown are the square root of the variance, /(l — f)/N, where 
/ is the radio- loud fraction and N is the number of quasars in the bin considered (e.g., Schneider 
et al. 1992). The impression from Fig. [6| is that for a given redshift bin, the radio-loud fraction 
is independent of optical luminosity. This result stands in contrast to the analysis of Goldschmidt 
et al. (1999) who find that the radio- loud fraction increases with luminosity for each redshift bin 
considered. Consideration of the radio-loud fraction plotted in Fig. 7 of their paper suggests that 
the optical luminosity dependence claimed at 1.3 < z < 2.5 depends largely on the poorly measured 
radio-loud fraction at Mb = —28. However, at 0.3 < z < 1.3, their data convincingly shows optical 
luminosity dependence in the radio-loud fraction. Comparing the ordinate between the two panels 
of Fig. |6] suggests that the radio-loud fraction remains approximately constant with redshift; we 
consider this next. 
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Fig. 7. — Fraction of radio-loud quasars as function of redshift. Two redshift ranges are plotted. 
No redshift dependence in the radio-loud fraction is evident. 



5.2. Radio-Loud Fraction as a Function of Redshift 



Since there is little evidence of optical luminosity dependence of the radio-loud fraction in our 
data set, we decrease our errors by considering the radio-loud fraction in a larger absolute magnitude 
range, —26 < Mb < —28. Fig. [7| shows the results of this analysis, with error bars calculated as 
above. At 1.75 < z < 2.5, we classify 20 out of the 153 quasars in the optical luminosity range 
considered as radio-loud, corresponding to 13.1 ± 2.7% of the quasars. At z > 4, we classify 4 out 
of the 34 optically-selected quasars as radio-loud (considering only the Schneider et al. (1992) and 
our targeted VLA survey). This corresponds to 11.8 ± 5.5% of the quasars being radio-loud. No 
evolution in the radio- loud fraction is detected. 

The FIRST/NVSS detections of optically-selected quasars at z > 4 also provides a lower limit 
to the radio-loud fraction at early cosmic epoch. Of the 107 z > 4 optically-selected quasars, 79 
have optical luminosities —26 > Mb > —28. From this restricted sample, 35 overlap with FIRST of 
which one was detected, implying a statistically unrobust radio-loud fraction > 3%. All 79 quasars 
in the luminosity range considered overlap with NVSS; four were detected, implying a radio-loud 
fraction > 5% at z > 4. 
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6. Conclusions 

We report on two programs to study the radio-properties of optically-selected quasars at high- 
redshift. First, we consider deep, targeted 5 GHz imaging of 32 z > 4 quasars selected from the 
Palomar multicolor quasar survey. Four sources are detected. We also correlate a comprehensive 
list of 134 z > 4 quasars, entailing all such sources we are aware of as of mid-1999, with two deep 
1.4 GHz VLA sky surveys. We find five new radio-loud quasars, not including one quasar identified 
in the targeted program. In total, we report on eight new radio-loud quasars at z > 4; only seven 
such sources are in the literature currently. 

We use this new census to probe the evolution of the radio-loud fraction with redshift. We find 
that, for — 25 > Mb ^ — 28 and 2 < z < 5, radio- loud fraction is independent of optical luminosity. 
We also find no evidence for radio-loud fraction depending on redshift. If the conventional wisdom 
that radio-loud AGN are preferentially identified with early-type galaxies remains robust at high 
redshift, this result could have implications regarding the formation epoch of late-type versus 
early-type galaxies. In hierarchical models of galaxy formation, one expects the late-type (less 
massive) systems to form first. Mergers are required to form the early-type (more massive) systems. 
Eventually, at high enough redshift, one would then expect the radio-loud fraction of AGN to fall 
precipitously. Our results show this epoch lies beyond z ~ 4, providing further evidence for an 
early formation epoch for early-type galaxies. 

We acknowledge the efforts of the POSS-II and DPOSS teams, which produced the PSS sample 
used in this work, and, in particular, N. Weir, J. Kennefick, R. Gal, S. Odewahn, R. Brunner, V. 
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